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1.  INTRODUCTION 

The  exact  solution  of  a  time-optimal  problem  for  a 
linear,  normal,  nt^-order  system  involves  an  (n-l)-dimen- 
sional  surface  as  a  switching  criterion  for  the  control 
variable.  Even  for  a  third-order  system  an  engineering 
realization  of  the  time-optimal  controller  is  very  com¬ 
plicated.  An  example  of  such  a  design  is  given  by  Athans 
and  Falb  (l).  For  more  practical  engineering  applications 
simpler  designs  for  higher  order  systems  were  developed 
using  suboptimal  controllers.  In  most  cases  the  great 
simplification  obtained  for  suboptimal  systems  well  com¬ 
pensated  for  the  small  deviation  from  optimality. 

The  basic  idea  used  in  most  suboptimal  designs  is 
that  the  higher  order  system  can  be  approximated  by  a 
second-order  system.  Among  others,  Kalman  used  this  ap¬ 
proach  and  derived  a  switching  curve  for  a  third-order 
system  (4).  His  method  essentially  consists  of  isolating 
the  system  modes  by  a  coordinate  transformation  and  then 
letting  the  two  dominant  roots  approximately  characterize 
the  dr  *ic  behavior  of  the  system.  Thus,  Kalman  implicitly 
a&  he  higher  order  system  has  two  roots  which 

are  mu<..  „  important  than  the  rest  of  the  roots  in  de¬ 

termining  the  behavior  of  the  system.  This  assumption  is 
the  starting  point  for  the  method  used  in  this  thesis  to 
analyze  higher  order  systems. 
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The  system  under  analysis  contains  small  parameters 
which  will  reduce  the  order  of  the  system  to  two,  when  they 
go  to  zero.  For  the  degenerate  second-order  system  sensi¬ 
tivity  functions  with  respect  to  these  small  parameters 
can  be  written.  Employing  these  sensitivity  functions,  a 
nearly  time-optimal  solution  is  found  by  deriving  a  simpli¬ 
fied  switching  function. 

The  object  of  this  thesis  is  to  find  a  method  to 
obtain  a  simplified  switching  function  for  a  higher  order, 
time-optimal  system  which  has  the  properties  described  in 
the  last  paragraph.  At  the  same  time  it  will  be  shown  how 
sensitivity  functions  can  be  used  to  find  the  approximate 
time  response  analytically. 
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2.  ANALYSIS  AND  DESIGN  PROCEDURE 

2.1  Formulation  of  Problem 

Consider  the  following  system  of  differential  equa¬ 
tions 

*1  =  x2 

x 2  «  klxl  +  ^2X2  +  x3 

^*3  =  -  x3  +  x^ 

*m*n  +  "  ^  +  u  C2"1) 

where  x^  are  state  variables,  and  k2  are  constants, 
are  small  parameters  and  u  is  the  control  variable. 
When  A j  =  0  for  j  =  l,2,--»,m  the  system  (2-1)  re¬ 
duces  to  the  second-order,  degenerate  system 

*1  “  x2 

X£  =  k^Xj,  +  k2X0  +  u  (2-2) 

Any  linear,  second-order  system  with  constant  coefficients 
can  be  represented  by  equations  (2-2). 


4 


For  the  second-order  system  consider  the  following 
problem.  Given  some  initial  conditions  and  the  constraint 
u  <1,  find  the  control,  u(t),  which  will  transfer  the 
state  of  the  system  to  the  origin  of  the  phase  plane  in 
minimum  time.  The  well-known  solution  to  this  problem  in¬ 
volves  a  switching  curve  that  divides  the  phase  piano  into 
two  halves.  Let  this  switching  curve  be  given  by 

fo(x1,x2)  =  0  (2-3) 

The  exact  solution  of  the  same  problem  for  the  n  *- 
order  system  (2-l)  involves,  as  mentioned  earlier,  (n-l)- 
dimensional  switching  surfaces  in  the  n-dimensional  phase 
space.  However,  it  is  reasonable  to  think  that  for  small 
>\j,  the  solution  of  the  n^-order  system  should  be  similar 
to  that  of  the  second-order  system.  In  other  words  for  a 
nearly  time-optimal  solution  for  the  n  -order  system  the 
switching  functions  can  be  expressed  in  terms  of  a  switching 
curve  in  the  x^x2-plane.  The  problem  is  to  find  this 
switching  curve. 


2,2  Sensitivity  Func t ions 


In  general  the  values  of  ^  could  be  different  as  long 
as  they  are  all  small  enough  to  make  the  exponential  tran¬ 
sients  negligible.  However,  there  is  not  much  loss  in 
generality  by  making  the  ^j's  all  equal  since  they  are 
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required  to  be  small  in  the  first  place.  What  •’small’*  means 
quantitatively  will  become  clear  later. 

The  state  variables,  x^(t,  7\.  )  for  i  »  1,2,  can  be 
approximated  by 


x^t,  7\  ) 


x^t.O)  + 


Sxt(t,  ) 

T* 


A=0 


(2-4) 


here  the  subscript  of  has  been  dropped  because  they  are 
now  assumed  to  be  all  equal.  Let  us  define  the  sensitivity 
functions  by 


wi  s 


dx?(t,7*  ) 


da 


A=0 


The  sensitivity  functions  are  obtained  from  the  following 
system  of  equations 


+  iLu  (2-5) 

b?\ 

There  are  two  difficulties  connected  with  the  solution 
of  these  equations.  First  of  all  the  function  w^(t)  is  not 
defined  at  t  =  0  in  all  cases  but  it  is  always  defined  for 
t  >  0,  For  this  reason  special  initial  conditions  have  to 
be  found  at  t  =  0+  where  the  "+"  sign  indicates  that  t  >  0 
by  a  small  amount.  This  problem  is  dealt  with  in  chapter  3. 


w 


2  »  klwl  +  k2w2 
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The  second  difficulty  occurs  at  the  switching  time, 
t8,  where  the  term  ^u/c)7i  is  discontinuous.  Consequently 
there  will  be  a  ''jump11  in  W2(t)  at  t  «  t8.  Let  us  denote 
this  "jump”  by  Aw2.  It  will  also  be  shown  in  chapter  3 
how 

Aw2  =  w2(ts+)  -  w2(ts")  (2-6) 

can  be  computed. 

Now  it  is  possible  to  write  the  sensitivity  functions 
for  all  t  >  0  and  t  4  ta#  Using  equation  (2-4),  approxi¬ 
mate  trajectories  of  the  n^-order  system  can  be  found. 


2.3  General  Form  of  Switching  Function 

Assume  that  the  switching  curve  for  the  n  -order 
system  is  defined  by 


f(x1,x2,?i)  «  0 


This  function  is  unknown  and  has  to  be  found  ip  some  way. 
An  approximate  expression  of  f(x^,x2,  )  in  terras  of  the 

switching  function  for  the  second-order  system  is  given  by 


f(x1,x2,  ?l  )  =  f(x1,x2,0)  + 


df  (x1}x2,  >|  ) 

b 


* 

7l»0 


(2-7) 


In  order  to  simplify  the  potation,  let  us  write 
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For  the  systems  considered  in  this  thesis,  it  is 
possible  to  express  the  sensitivity  functions  as  a  linear 
function  of  the  x^*s  of  the  second-order  systems.  In 
general 

wj^  =  alxi  +  a2x2  +  a3 

w2  =  blxl  +  ^2X2  +  ^3  (2-11) 

where  the  a's  and  b's  are  constants.  These  constants  de¬ 
pend  on  initial  conditions  in  most  cases.  Since  a  switching 
curve  should  be  independent  of  initial  conditions,  it  is 
not  meaningful  to  use  the  constants  of  (2-11)  in  the 
switching  function.  The  usefulness  of  equations  (2-11)  lies 
in  the  fact  that  they  give  the  correct  functional  depend- 
ence  on  the  state  variables,  x^,  for  the  approximate 
switching  curve  (2-9)  regardless  of  the  initial  conditions. 

From  (2-9)  and  (2-11)  the  switching  function  for  the 
ntll-order  system  is  thus  determined  vithin  some  undefined 
constants.  These  unknown  constants  will  be  computed  by  some 
other  means.  But  first  let  us  consider  the  problem  of 
finding  the  special  initial  conditions  for  the  sensitivity 
functions. 
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3.  SPECIAL  INITIAL  CONDITIONS  FOR  SENSITIVITY  FUNCTIONS 
3„1  Changes  Method 

It  was  mentioned  that  a  difficulty  arises  in  finding 
the  initial  conditions  for  the  sensitivity  equations  (2-5),, 
The  function  w^(t)  may  not  be  continuous  at  t  *  0o  To 
illustrate  this  point  consider,  the  simple  system 


=  -X2  +  u  (3-1) 

with  zero  initial  conditions  and  u  as  a  step  function. 

When  7\  a  0,  the  sensitivity  equation  is 

wx  =  0  (3-2) 

The  usual  initial  conditions  for  sensitivity  functions  with 
respect  to  variable  parameters  are  zero.  But  this  will  make 
w^(t)  =  0  which  is  obviously  wrong.  The  trouble  is  that 
(3-2)  is  not  defined  at  t  =  0.  One  way  to  avoid  this  diffi¬ 
culty  is  to  use  w^(0‘*’)  as  an  initial  condition  for  (3-2), 

A  direct  solution  of  (3-1)  gives 

xL  »  ut  +  u2i(e”t/?'  -1)  (3-3) 

Differentiating  (3-3)  with  respect  to  and  then  setting 
?i  s  0,  we  have 
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=  w^(t)  =  -u  (3-4) 

=0 

Therefore  *  correct  initial  condition  for  (3-2)  is 
w1(0t)  o  -  u, 

Chang  developed  a  method  to  find  the  initial  condi¬ 
tions,  w^(0+),  (2).  For  the  system  (2-1)  this  method  applies 
in  the  following  way.  First  the  Laplace  transforms  are 
written  for  all  equations  in  (2-1),,  All  initial  conditions 
are  equal  to  zero  except  for  x^(0)  and  X2(0).  But  x^(0)  and 
x2(0)  can  also  be  set  to  zero  because  they  do  not  affect 
the  end  result.  Solving  the  transform  equations  for  Xi(b,/\  ) 
and  X2(s,  )  gives 

X^(s,  )  =  _ u(0) 

s(s2-k2s-k^)  (?1^  s+l)---(  ms+l) 

(3-5a) 


X2(^s^)=  _ njO) _ _ 

(s2-k2s-k1)  (  ?i  j^s+l)-— (  7\  ms+l) 

(3-5b) 


Remember  that  }\  =  ^  for  j  =  l,2,---,m0  Differentiating 

(3-5a)  with  respect  to  ^  and  then  letting  ^  =  0,  we  get 


C>X1(s,  >  ) 


=  WL(s) 

*=0 


-mu(0) 

1,,M-  - 

s  -kns-k 


1 
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The  initial-value  theorem  gives 

w,(0+)  =  lim  sW.  ( 8 ) 

A  s-^£x>  L 

=2  0 

Now  differentiating  (3-5b)  with  respect  to  ^  and  setting 
^  a  0  gives 

=  Wjj(s)  =  -mu(0)s 

2 

7l=0  s  -k2S-k^ 

Again  from  the  initial-value  theorem 
wo(0*)  sa  lim  sW«(s) 

Z  S**-C*3  *■ 

-  -mu(0). 

In  summary  the  initial  conditions  for  the  sensitivity 
functions  defined  by  the  equations  (2-5)  are 

Wl(0+)  =  0 

w2(0+)  =  -mu(0)  (3-6) 

3.2  Kokotovlc^Rutman  Structural  Method 

An  interesting  alternative  approach  to  determine  the 
initial  conditions  of  the  sensitivity  functions  is  given  by 
Kokotovic*  and  Rutman  (5),  The  method  consists  of  drawing 
a  block  diagram  of  the  original  system  and  the  sensitivity 
system  for  ^  >  0  and  then  manipulating  the  blocks  until  it 
is  possible  to  let  =  0. 


^X2(s,  >  ) 
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To  illustrate  this  method,  let  us  find  the  initial 
conditions  of  the  sensitivity  functions  for  the  following 
system  which  will  be  used  in  an  example  later. 


*1  =  x2 


*2  =  x3 


£3  a  -X3+U 


XL(0)  a 


x2(0)  =  r2 
x3(0)  =  0  (3-7) 


The  corresponding  system  of  sensitivity  equations  is  of 
the  form 


W1  =  W2 


*2  =  W3 


w3  =  -w8 


-  Xo  + 


Bu 


(3-8) 


In  Figure  la  the  original  system  diagram  and  the  sensiti¬ 
vity  model  are  shown.  Figure  lb  is  obtained  after  two 
transformat'.  >ns.  First  the  point  of  application  for  x3  is 
moved  to  the  other  side  of  the  integrator  and  a  differential 
operator  is  added  in  the  connecting  branch  of  the  two 
models.  Now  it  is  possible  to  let  7\  =  0.  The  transfer 
function  from  u  to  x3  as  ?\  goes  to  zero  is 


lim  1 

^s+1 


=  1 
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The  final  diagram,  Figure  lc  has  the  differential  operator 
removed  by  moving  the  point  of  application  to  the  other  aide 
of  the  integrator  in  the  sensitivity  model. 

From  Figure  lc  the  initial  conditions  can  be  written 
by  inspection. 

wl(0+)  =  0 

W2(0+)  ss  -u(0) 

3.3  Discontinuity  at  Switching  Time 

The  function  ^(t)  is  discontinuous  at  t=t8  because 
the  control  variable,  u(t),  is  switching  from  one  limiting 
value  to  the  other  one  instantaneously.  As  mentioned  pre¬ 
viously  there  will  be  a  ''jump”  in  W2(t)  at  t»ts  and  it 
becomes  necessary  to  find 

w2(ts+)  a  w2(tg+)  +  Aw2  (3-9) 

in  order  to  solve  the  sensitivity  equations. 

If  it  is  possible  to  find  a  transformation  of  the 
system  (2-1)  su^h  that  the  transformed  system  is  a  func¬ 
tion  of  t-t8  but  otherwise  has  the  tame  form  and  initial 
conditions  as  the  original  system,  the  condition  (3-9) 
can  be  determined  from  the  initial  conditions  of  the  trans¬ 
formed  system.  First  let  us  show  that  such  a  transformation 
exists  for  the  system  (2-1), 
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The  system  (2-1)  can  be  expressed  in  matrix  form  as 


shown  below. 


x(t)  =  A  x(t)  +  z 


(3-10) 


when  x  and  z  are  n-vectors  and  A  is  an  (nxn)  matrix.  The 
matrix  A  is  given  by 


A  = 


:1  k2  1  Q 

0  -L  1 

X*  N  X 

\  N  \ 

\  N  \ 

~  o  -I  1 

0  » 

0  -i 

a 


(3-11) 


The  vector  z  has  only  one  non- zero  element  in  the  n^  row 
which  is  u/?i  . 

Let  us  define  the  new  time  variable,  T,  as 


T  »  t  -  t„ 


(3-12) 


The  initial  conditions  for  (3-10)  are  assumed  to  be  zero. 
Suppose  now  that  there  exists  the  following  matrix  equa¬ 
tion 


i(T)  =  A  y(T)  +  z 


(3-13) 


with  the  initial  conditions  y(0)  «=  0.  The  next  step  is  to 
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find  the  transformation  of  x  which  will  give  2  (3-13). 

The  general  solutions  of  (3-10)  and  (3-13)  are 
,t 

A(t-p) 

(3-14) 

0 

"T  A(T-d) 

(3-15) 

0 


A(t-p) 

x(t)  =  /  e  z  dp 


,  %  .  A(T-p) 

i(T)  a  /  e  z  dp 


In  the  above  equations  p  is  a  dummy  variable.  Substituting 

t  a  t,  +  T  into  (3-14)  gives 
8 


,T+t. 


A(t  +T-p) 


x(T)  =  /  e  z  dp 

<3 

,T+t 


At 


,  3  A(T-p) 

=  e  J£(T)  +  I  e  z  dp  (3-16) 

T 


Solving  (3-16)  for  ^(T),  we  have 


Z(T) 


-At 

e  sx(t)  + 


A(T-p) 
e  z  dp 


T+ts 


(3-17) 


Equation  (3-17)  gives  the  transformation  from  x(t)  to 
^(T),  so  that  ;jr(T)  satisfies  the  equation  (3-13),  As  long 
as  the  state  transition  matrix  exists,  it  is  possible  to 
find  this  transformation. 
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relation  with  respect  to  7l  %  we  get  the  equations  relating 
Vj^T)  and  w^t). 

Consider  the  particular  case  when  k^  **  kg  °  0*  The 
solution  of  equation  (3-16)  for  this  case  is 

ut,2 

Xi<t)  a  y^(T)  +  tsy2(T)  +  -7— 

X2(t)  a  y2^T)  ♦  uts  (3-21) 

In  equations  (3-21)  u  is  considered  to  be  constant  and 
continuous.  The  fact  that  it  is  actually  discontinuous  at 
the  switching  time  is  accounted  for  by  the  special  initial 
conditions  of  the  transformed  system  (3-13),  With  this  in 
mind  (3-21)  can  be  differentiated  with  respect  to  7\  and 
the  result  is 

w1(t)  =  V1(T)  +  t8v2(T) 

w2(t)  «  V2(T)  (3-22) 

Evaluating  (3-22)  at  t  =  t8+  or  T  =  0+  and  substituting 
(3-18)  gives 

wi(t8+)  =  -  tgmu(ts+) 

W2(ts+)  a  -  mu(tg+) 

Since  u(t5+)  =  -u(t=0)  equations  (3-23) become 
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Wj.Ctg*)  =  t8mu(t=0) 

W2(tg+)  =  mu(t=0)  (3-24) 

For  this  particular  system  the  “jump"  in  at  the 
switching  time  is 

A  \*2  =  -  2mu(t=0)  (3-25) 

In  a  similar  way  the  "jump"  conditions  can  be  derived  for 
other  cases  of  the  general  system  (2-1 ). 


\ 


4.  SWITCHING  FUNCTION  FOR  A  PARTICULAR  SYSTEM 
4.1  Functional  Form 

The  general  form  of  che  system  under  consideration  is 
given  by  (2-1).  When  =  k2  =  0,  we  have  the  following 


system 


fcl  = 


x2  =  x3 


*  -xw  :■  j. 


^  *n  =  -xn  +  u 


(4-1) 


The  corresponding  second-order  system  for  *  = 


0  is 


*1  =  x2 


c2  =  u 


(4-2) 


Let  the  initial  conditions  be  denoted  by 
xl(0)  =  r^ 

x2(°)  =  r2 

The  solution  of  (4-2)  is 


(4-3) 
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(4-8) 
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Differentiating  (4-8)  with  respect  to  x^  and  gives 

ltn 


ax. 


=  1 


l£c 

2>X<- 


x. 


(4-9) 


If  we  substitute  (4-7)  and  (4-9)  into  (2-10)  we  get 


fl  =  (a2x2  +  a3^  ^  +  |x2 I  b3  ^  (4-10) 

Defining  new  constants  (4-10)  can  be  written  as 


i1  =  b  x2  +  c> 


where 


b 


a2 


(4-11) 


c  =  a3 

Finally  the  total  switching  function  from  (4-8)  and  (4-11) 
is 

f  =  xl  +  I  x2  | x2  |  +  b^x2  +  (4-12) 

The  next  task  is  to  determine  the  two  constants  b  and 
c  uniquely  and  independent  of  the  initial  conditions.  This 
problem  is  solved  in  the  next  section. 

The  procedure  outlined  above  can  easily  be  applied  to 
two  other  cases  of  the  general  system  (2-1):  1.)  when 
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k^  +  0  and  kg  4  0  2.)  k^  4  0  and  kg  4  0  and  the  eigen¬ 

values  of  (2-2)  are  distinct  and  real.  For  the  case  of 
complex  eigenvalues  difficulties  arise  because  the  optimal 
switching  curve  cannot  be  expressed  in  a  nice  analytical 
form.  However,  if  the  optimal  switching  curve  is  approxi¬ 
mated  by  some  analytical  function  of  the  state  variables, 
it  is  possible  to  treat  this  case  in  the  same  way. 

4.2  Determination  of  Constants 

In  order  to  compute  the  constants  b  and  c,  it  is 
necessary  to  derive  the  equations  for  the  approximate 
trajectory  of  the  nth-order  system  as  given  by  (2-4).  In 
this  derivation  the  following  initial  condition  will  be 
used. 


x1(0)  =  rl 

Xg(0)  =  0  (4-13) 

The  initial  conditions  (4-3)  could  also  be  used.  Since  the 
final  result  is  independent  of  the  initial  conditions,  it 
is  more  convenient  to  use  (4-13).  The  mathematical  analysis 
is  simpler  with  (4-13).  However,  the  initial  conditions 
have  to  be  general  enough,  so  that  every  point  on  the 
switching  curve  can  be  reached  through  them.  This  require¬ 
ment  is  satisfied  by  (4-13). 

Another  simplification  of  the  analysis  is  obtained  by 
only  considering  one  half  of  the  switching  curve.  The  other 
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half  Is  determined  by  sy^etry.  By  restricting  r,  to 
poaitivo  values,  the  trajectories  will  only  reach  che  lower 
half  of  the  switching  curve. 

For  this  case  the  optimal  control  sequence  for  the 
system  (4.2)  is  u  (.2,  The  solutlon 


For  0  <  t  <  t  • 
-  -  us* 


xi  -  t  L 

1  “  T  +  ri 


x2  =  -t 
For  t  <  t  < 

S  —  —  • 


(4-14) 


xi  =  r  •  2v  *  *  rx 


x2  ■  t  -  2t, 


(4-15) 


The  solution  of  the  corresponding  sensitivity 

tions  (4.5)  with  the  initial  conditions  (3-6)  and  the 
"jump"  condition  (3-25)  is 


equa- 


For  0  <  t  <  t  : 


W1  =  mt 


Wr 


m 


(4-16) 


For  t8  <  t  <  t£ : 
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w^  a  -  mt  +  2mtQ 

W2  =  -  m  (4—1 7 ) 

t*u 

The  approximate  equations  for  the  nc  -order  system 
(4-1)  are 


For  0  <  t  <  tc: 

.  t2 

x,  =  +  m/it  +  r, 

1  2  1 

x2  a  -t  +  m7\  (4-18) 

For  t8  <  t  <  tf : 

t2  0 

x^  =  y-  -  (2tg  +  m/\)t  +  ts  +  2m>its  +  r^ 

X£  +  t  -  2ts  -  m  A  (4-19) 


An  expression  for  the  switching  time,  tg,  is  obtained 
from  the  intersection  of  the  switching  curve  (4-12)  and 
the  trajectory  defined  by  equations  (4-18).  The  final 
result  is 


t 


(b-2)m>\ 


+ 


+ 


2j*(b2+2)m2  ?\2  +  cm?l 


(4-20) 


A  condition  that  has  to  be  satisfied  by  the  trajec¬ 
tory  (4-19)  is  that  it  should  reach  the  origin  at  the 
final  time,  tf.  In  other  words,  x^(tf)  =  x2(tf)  =  From 
equations  (4-19)  this  will  only  happen  if 
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=  2tg  + 


(4-21) 


and 


12  2 

-  -  +  t_  +  2m  +  r<i  =  0 

9  r  s  si. 


(4-22) 


Substituting  (4-21)  into  (4-22)  and  solving  for  t8  gives 

% 

(4-23) 


fcs  = 


e  _m2*2 

rl  —2~ 


The  two  expressions  for  ts,  (4-20)  and  (4-23),  are 
equal  if 


b  =  2 

c  =  -2m ^ 


(4-24) 


It  turns  out  that  this  is  the  only  choice  of  the  constants 
b  and  c  which  is  independent  of  r^  and  at  the  same  time 
makes  equation  (4-20)  equal  to  (4-23), 

Putting  the  constants  (4-24)  into  the  switching  func¬ 
tion  (4-12)  and  correcting  for  the  sign  gives 


f  = 


x. 


,  1 

+  1  x2 


+  2m  ^  X£  +  2m  ^ 


(4-25) 


Remember  that  m  =  n-2.  So  for  a  third-order  system  m  =  1. 
The  switching  function  for  that  case  is  almost  in  perfect 
agreement  with  one  derived  by  Kalman  for  the  same  system 
by  a  completely  different  method  (4),  The  only  difference 
is  that  in  (4-25)  the  constant  multiplier  of  ^  is  2  while 
lyi  Kalmans  equation  it  is  1/2, 
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Figure  2  shows  two  approximate  switching  curves  for  a 
third-and-fourth-order  system  when  >  =  1/2.  The  optimal 
switching  curve  for  the  second-order  system  is  also  shown 
for  comparison.  Note  that  the  two  approximate  switching 
curves  do  not  go  exactly  through  the  origin.  The  "miss” 
at  the  origin  is  proportional  to  7\  If  7\  is  kept 
smaller  than  one,  the  curves  will  come  quite  close  to  the 
origin. 


In  this  section  the  time  response  of  the  third-order 
system  (3-7)  will  be  computed  for  two  different  initial 
conditions  and  7\  a  0,5  using  the  derived  switching  func¬ 
tion  (4-25), 

First  let  us  choose  r^  =  4,125  and  r2  -  0.  Conse- 
wuently  the  switching  time,  tg,  is  2  and  the  final  time, 
tf,  is  4.5,  The  equations  describing  the  trajectory  are 
given  below. 


For  0  <  t  <  ts: 


w2  =  1 

x^  a  -  0.5t^  +  4.125 


\ 
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xj_  +  0.5w^  =  -  0.5t^  +  0.5t  +  4.125 
+  0.5w2  =  -  t  +  0.5 
For  tg  <  t  <  tf: 

w^  «*  -  t  +  4 
w2  a  -  1 

x^  ts  0.5t^  -  4t  +  8.125 
x2  =  t  -  4 

x^  +  O.Sw^  a  0.5t^  -  4.5t  +  10.125 
x^  +  0.5w2  =  t  -  4.5 

Let  the  exact  solution  of  the  system  of  differential 
equations  with  u  switching  form  -1  to  +1  at  t  =  tg  be 

•ff 

denoted  by  x^  for  i  =  1,2.  Figure  3  shows  the  various 
curves  for  this  example.  Note  that  there  is  an  error  of 
about  ?[  ^  between  the  exact  and  approximate  curve  of  x^ 
at  t  s  tf.  Since  the  optimal  time  for  the  second-order 
system  when  =  0  is  4.06,  one  can  see  that  the  time 
response  is  very  nearly  time-optimal. 

As  a  second  example  let  us  use  the  initial  conditions 
r^  «  0  and  =  2,  The  constants  for  the  switching  curve 
were  derived  with  zero  initial  conditions  for  X2.  This 
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example  will  illustrate  the  fact  that  the  switching  curve 
nevertheless  is  valid  for  non-zero  initial  conditions  of 
x2* 

The  formulas  for  tfl  and  tf  are  easily  derived  in  the 
same  way  as  before.  The  result  is 

ts  a  1.707r2  +  0.707 

tf  a  2.414  (r2  +  ^  ) 

For  this  example  ts  «  3.768  and  tf  =  6.035. 

The  trajectory  is  described  by  the  following  equations. 
For  0  <  t  <  ts: 

W]_  =  t 

w2  a  1 

xl  =  -  0.5t2  +  2t 
X2  «  -  t  +  2 

xL  +  0.5Wl  =  -  0,5t2  +  2.5t 
x2  +  0.5w2  a  -  t  +  2.5 
For  t8  <  t  <  tf : 

=  -t  +  7.536 
w2  =  -  1 
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x^  a  0.5t^  -  5.535t  +  14.20 
X2  «  t  -  5.535 

Xj^  +  0.5*^  a  0.5t^  -  6,035t  +  17.97 
X2  +  0.5w2  -  t  -  6.035 

Figure  4  shows  the  various  curves  for  this  example.  Again 
the  time  response  is  a  good  approximation  of  the  time- 
optimal  solution. 
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5.  SUMMARY  AND  CONCLUSION. 

A  method  has  been  presented  to  derive  simplified 
switching  functions  for  the  design  of  third  or  higher  or¬ 
der,  time-optimal  control  systems.  They  type  of  system, 
that  has  been  considered,  is  characterized  by  small  para¬ 
meters  which  increase  the  order  of  the  system.  The  sensi¬ 
tivity  functions  with  respect  to  these  small  parameters 
have  been  the  main  tool  of  analysis.  Two  existing  methods 
have  been  used  to  determine  special  initial  conditions  for 
the  sensitivity  equations  because  these  equations  are  dis¬ 
continuous  at  t  »  0.  It  has  been  shown  that  the  same 
methods  can  be  applied  to  deal  with  the  discontinuity  of 
the  sensitivity  functions  at  the  switching  time  by  a  trans¬ 
formation  of  the  original  system. 

A  particular  system  has  been  analyzed  in  more  detail. 
The  degenerate,  second-order  system  for  this  case  is  the 
well-known  double  integral  plant.  The  complete  switching 
curve  has  been  determined  for  this  system.  Using  this 
switching  curve  the  approximate  time  response  was  computed 
with  the  aid  of  sensitivity  functions.  The  results  show 
that  the  response  is  nearly  time-optimal. 

For  a  system  characterized  by  a  second-order  system, 
which  is  more  complicated  than  the  double  integral  plant, 
it  might  be  difficult  to  analytically  determine  unique  con¬ 
stants  in  the  switching  function.  <In  that  case  perhaps  the 
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best  constants  can  be  found  experimentally.  As  already  in¬ 
dicated,  complex  roots  of  the  second-order  system  also  pre¬ 
sent  special  problems.  All  these  cases  need  further  study, 
for  which  there  was  no  time  in  connection  with  this  thesis 
unfortunately. 

This  investigation  above  all  demonstrates  how  sensi¬ 
tivity  functions  with  respect  to  small  parameters,  that 
change  the  system  order,  can  be  used  to  greatly  simplify 
the  analysis  of  higher  order  systems. 
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